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1. Introduction 

There are current experimental and theoretical activities to investigate the decays of heavy- 
light hadrons. Study of their weak decays is of interest for determining the CKM matrix of quark 
mixing, and there are theoretical and experimental efforts to test and overconstrain the Standard 
Model and find signatures of New Physics. 

It is of interest to study charmed hadrons on fine lattices where discretization effects are very 
small. While unquenched data from very fine lattices is not yet easy to obtain, the QCDSF col- 
laboration has undertaken a quenched calculation of charmed and also bottom decay constants, 
weak matrix elements and meson spectra on a quenched lattice with an inverse lattice spacing of 
a -1 ~ 5 GeV. Results on heavy-light and light meson decay constants have been presented in a 
previous paper [jlj]. A particularly interesting result of this calculation was a relatively low value 
of the decay constant fo s . We have analyzed further matrix elements and would like to present 
our findings here. In Section |2| we give a short description of simulation and analysis details. Our 
calculation of the pseudoscalar decay constants using two different interpolating operators for the 
pseudoscalar meson is discussed in Section [3} In Section || we give preliminary results for form 
factors of semileptonic decays of pseudoscalar heavy-light mesons to pseudoscalar light mesons. 

2. Simulation parameters 

Our quenched lattices are generated using the Wilson gauge field action at j8 = 6.6. The 
lattice spacing, determined using the Sommer parameter ro = 0.5 fm from [0], is a = 0.040 fm or 
aT x = A. 91 GeV. The lattice size is 40 3 x 80. For the results reported here, we have analyzed 114 
gauge field configurations. 

We use 0(a) improved Wilson quarks. The value of the clover coefficient has been determined 
nonperturbatively by Ref. [Q]. We work with three 'light' and four 'heavy' hopping parameters. 
Their values and the corresponding quark masses can be found in Table [j]. 

K 0.13519 0.13498 0.13472 0.13000 0.12900 0.12100 0.11500 

am q 0.0076 0.013 0.020 0.14 0.16 0.31 0.37 

M[GeV] 0.55 0.69 0.87 2.8 3.2 5.6 7.2 



Table 1: Hopping parameters, the corresponding 0(a) improved quark masses (see Eq. (2.1)) and approxi- 
mate pseudoscalar meson masses (denoted as M). 

In our analysis we parameterize the light quark masses using the 0(a) improved quark mass 

am q = (l+b m am q ), (2.1) 

where am q = \(\ — The nonperturbatively determinated value of b m is taken from Ref. [Q]. 
As a measure for the heavy quark masses we use the mass of the heavy-light meson, where the 
light quark mass is extrapolated to the average of the u and d quark mass. 

We calculate matrix elements of the 0(a) improved axial vector current containing the quarks 
q\ and q 2 : 

A' 4 =Z A (\+ b A am q ) (A 4 + ac A d A P) , (2.2) 
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where A^(x) = q { (x) 74 75*72 (X) is the local axial vector lattice current, and P(x) = q l (x)y= > q2{x) 
the local pseudoscalar density. With we denote the symmetric lattice derivative. We use the 
nonperturbatively determined values for ca from ^ and for Za from [|[]. For the coefficient bA of 
the bare quark mass am q , we use the value of Ref. calculated in one-loop perturbation theory, 
and the boosted lattice coupling. 

The 0(a) improved vector current is given by 

VjJ =Z v (\+am q b v ) {Vfj, + iac v d v T^ v ) , (2.3) 

where = q x y^qi and T^ v = q 1 ^ [7^ , y v ]q2- We use for Zy and by the nonperturbatively calculated 
results from For the coiTection due to the tensor current we have only made a preliminary 
order-of-magnitude estimate, using a rational interpolation ^ of the nonperturbative values for cy 
from [g]. 

3. Pseudoscalar decay constants 

We calculate the decay constant / of pseudoscalar mesons M at zero momentum from 

f=L(0\A{\M). (3.1) 



Previously ([|I|]) we used operators of the type P to project at the source onto the pseudoscalar 
meson. Since different correlation functions may be subject to different systematic errors, we 
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Table 2: Pseudoscalar heavy-light meson masses and decay constants at O(a ) from (A4P) and (A4A4) 
correlation functions at the simulated hopping parameters. 

would like to determine the matrix elements of A4 also using the temporal component of the axial 
vector current operator, i.e. from correlation functions of the form 

ClWC) = X>ift t)Af (0)>. (3.2) 
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The index i stands for local (i = L) or Jacobi smeared (i = S) operators. If we write the decay 
constant as 

f = Z A (\+b A am q )(fW+ac A f^), (3.3) 

we find the unimproved contribution /W from the correlation function in Eq. (3.2). 

We extract masses and amplitudes from single state fits to the correlation functions. The 
pseudoscalar meson masses and the results for /(°) are compared for heavy-light mesons at the 
simulated quark masses in Table |2| with the corresponding values from The results agree 
within the statistical errors. 
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Table 3: Results on decay constants from different correlators. The first error is statistical, the second error 
on the results from [jl|] is systematic and the third from the uncertainty in the experimental value of ro. We 
estimate the systematic errors in the new results to be similar. 

We fit the matrix elements to linear functions in the light quark masses and, where applicable, 
to quadratic functions in the heavy quark masses, using the method described in [|l|]. To determine 
the physical light, strange and heavy quark masses we also use the same method as for the central 
values given in [jl|]. The decay constant results are listed in Table |3[ We take over the errors for the 
heavy-light decay constants from pi], including uncertainties from setting the quark masses to the 
physical values, discretization effects, errors in the nonperturbative renormalization constants and 
from finite volume effects. A 10% uncertainty in the experimental value of ro is also included. We 
did not yet perform the same reanalyses for the data from the A4A4 correlators and just assume that 
the systematic errors are very similar. We did not attempt an estimate of the systematic error for 
the light meson decay constants. 

Our result for fo s is still smaller than the recently updated experimental value of the D s decay 
constant of 275(10)^(5)^ MeV by CLEO [0. 

4. Semileptonic form factors 

We describe preliminary results for matrix elements of semileptonic decay of pseudoscalar 
heavy-light mesons ('D') with mass M/, and momentum p^, to pseudoscalar light mesons ('71') with 
mass Mi and momentum pi. The decay goes via a vector current (x) with the matrix element 

JK ll = {%{pi)\V lx {Q)\D{p h )). (4.1) 
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We would like to determine the two form factors which are often used to parameterize the matrix 
element 

= (p h + Pi -Aq)nF+(q 2 )+ Agrafe 2 ), (4.2) 
with q = Ph~ Pi and A = (M 2 — Mf)/q 2 Jl3|]. The matrix element can be extracted from the three 
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Figure 1: 'Divided' three point functions at the hopping parameter of the heavy quark ff/, = 0.12900 and 
the hopping parameter K = 0.1 3498 for the light and spectator quark. The meson momenta are /?/, = ( 1 , 0, 0) 
and q = (1,0,0) in units of 2n/L. The crosses denote the data points, the solid lines the central values and 
the dashed lines the error bars of the fit. On the left we show | ^#4 1 and on the right | J£ 1 1 . 

point correlation function of a pseudoscalar light meson at time zero, a vector current at time t x and 
a pseudoscalar heavy dight meson at time t y = T/2, where T is the time extent of the lattice: 

4 3) (t x ,t y ) = £ e-^V^(if (y, t y )V, (x, t x )Pf(0))- (4-3) 

Pf and Pf are Jacobi smeared operators with the structure qjsqi and q h ysq. With qt, we denote the 
heavy quark, with q\ the light quark, and with q the spectator quark. The vector current is of the 
type V M = qflyflh. 

The large time behavior (0 <C t x <C t y ) or (0 <C t y <C t x ) of the three point function is given by 

Cf{t x ,t y ) — > ^L^Le-^e-*^ (K( Pl )\V,\D(p h )),t x <t y , 
Cf\t x ,t y ) — > ±^Le-W-t*) e -m*-h) ( K ( Pl )\ V ,\D(p h )),t x >t y , 

where the prefactors are given by Zf = \(§\P h s \D{ph))\ and Zf = \(0\Pf\n(pi))\. The correlation 
functions are symmetric or antisymmetric around T/2. We bin over the contributions from t x > t y 
and t x < t y , and divide the three point functions by the prefactors extracted from fits to smeared- 
smeared two point functions in a bootstrap loop. The 'divided' three point function is fitted to a 
constant to extract the matrix element. An example is shown in Fig. [j] In this preliminary analysis 
we considered in particular momentum combinations where the D and the n both have spatial 
momenta < 2n/L, where L is the spatial lattice extent, and the momenta are aligned, i.e.: 

p h = (0,0,0), -q = (0,0,0), (1,0,0), (0,1,0), (0,0,1), (-1,0,0), 
p h = (1,0,0) , -q = (0,0,0) , (-1,0,0) , (-2,0,0) , 
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in units of 2%jL. The kinematic conditions are such that for each heavy and light meson mass, each 
different combination of pf, and q leads to a different value of q 2 = (£/, — £/) 2 — q 2 . To extract the 
form factors for each given value of q 2 separately, without making a guess for the q 2 dependence 
of the form factors, we only have the four Lorentz components of the matrix elements at our 
disposition. 

We have exploited the fact that the form factors appear in the spatial and temporal components 
of with different momentum prefactors, so we obtain a linear equation system which we can 
solve for Fo and F + . 

Results for the form factors at several heavy K values and with the light and spectator quark 
mass close to the strange quark mass are shown in Fig. ||. For a preliminary estimate of the statis- 
tical error we use error propagation of the bootstrap errors of the spatial and temporal components 



of jtf^. In the results shown in the Figure, the 0(a) correction to the vector current (see Eq. (|2.3[)) 
is not yet included. We have however estimated its magnitude to be at most a few percent for heavy 
quark masses around the charm. Although the quark masses are not tuned to the physical values 
matching K mesons, and the results correspond to decays of D s to pseudoscalar strange mesons, 
we attempt to compare them to experiment and other lattice results. For k = 0.12900, close to 
K:harm> we find Fq(0) « f+(0) si 0.75. For D —> K decays, a recent experimental determination 



from BABAR has found a value of 0.727(11) [11], a lattice calculation with three flavors of dy 



namical quarks on a coarser lattice using Fermilab heavy quarks quotes F + (0) = 0.73(8) [12], and 
a quenched calculation on a coarser lattice quotes F + (0) = 0.66(4) [O]. So there seems to be at 
least a rough agreement, but for a more precise result the tuning of quark masses to the physical 
values will have to be performed. 
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Figure 2: Form factors for semileptonic decays for various heavy quark masses. The light and spectator 
quark hopping parameters are fixed to the K value 0.13498. Filled symbols denote Fo, open symbols F+. 
Squares correspond to K h = 0.129 {M h — 2.1 GeV), circles to K h = 0.121 (M h — 3.4 GeV), and triangles to 
K h = 0.115 (M h =4.2 GeV). 
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